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Abstract 

Although  significant  advances  have  been  made  in  hy¬ 
personic  boundary-layer  transition  prediction  in  the 
last  several  decades,  most  design  work  still  relies  on 
empirical  correlations  or  wind  tunnel  tests.  Codes 
using  the  semi-empirical  method  will  need  to  be 
verified  and  validated  before  being  used  for  expen¬ 
sive  flight  vehicles.  The  STABL  code  package  and  its 
PSE-Chem  stability  solver  are  used  to  compute  first 
and  second  mode  instabilities  for  both  sharp  and 
blunt  cones  at  wind  tunnel  conditions,  with  laminar 
mean  flows  provided  by  the  DPLR2D  Navier-Stokes 
code.  Stetson’s  3.81  mm  blunt  cone  case,  a  sharp 
cone  at  Mach  3.5,  and  a  very  blunt  cone  at  Mach 
8  are  analyzed.  The  computed  transition  locations 
agree  well  with  previous  computations  by  other  re¬ 
searchers,  but  larger  differences  cire  seen  in  the  local 
ampliflcation  rates  for  the  Stetson  blunt  cone  case. 
The  applicability  of  various  transport  property  mod¬ 
els  cind  their  effect  on  boundary  layer  stability  are 
examined.  This  work  helps  to  extend  the  applicabil¬ 
ity  of  STABL  to  low-temperature  flows. 


Nomenclature 


C„ 

Specific  heat  at  constant  volume 

E 

disturbance  kinetic  energy 

Gv 

vibrational  energy 

k 

thermal  conductivity 

M 

Mach  number;  molecular  weight 

P 

pressure 

Reynolds  number  Re  = 

Tn 

nose  radius 

S 

arc  length  along  the  body  surface 

T 

temperature 
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t 

time 

u,v,w 

velocity  in  x,y,z  directions 

X 

distance  in  streamwise  direction 

y 

distance  in  body-normal  direction 

z 

distance  in  spanwise  direction 

a 

streamwise  complex  wavenumber 

P 

spanwise  complex  wavenumber 

first  coefficient  of  viscosity 

V 

kinematic  viscosity 

P 

density 

<T 

disturbance  growth  rate 

4> 

vector  of  disturbance  quantities 

streamwise  computational  coordinate 

shape  function;  wave  angle 

n 

domain  of  integration  for  PSE 

CJ 

frequency 

Superscripts 

mean  component 
fluctuating  component 

Subscripts 

i  imaginary;  body-parallel  direction 

j  body-normal  direction 

o  stagnation;  critical  location 

rot  component  of  the  rotational  mode 

s  specific  chemical  species 

tr  translational;  transition 

vib  vibrational 

00  fireestream  conditions 

Introduction 

Despite  more  than  fifty  years  of  hypersonics  re¬ 
search,  many  critical  areas  are  still  poorly  under¬ 
stood  [1].  Hypersonic  flight  remains  a  top  priority 
in  defense  research,  as  any  vehicle  that  ascends  to 
or  firom  space  must  traverse  the  hypersonic  regime. 
The  reduction  in  forward  deployed  forces  since  the 
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end  of  the  Cold  War  has  increased  the  need  for  hy- 
personics,  as  a  need  remains  to  be  able  to  strike  high- 
value,  time-critical  targets  anywhere  in  the  world  in 
minutes  or  hours  [2].  Hypersonic  flight  is  a  key  com¬ 
ponent  of  most  plans  to  provide  that  capability. 

The  location  and  extent  of  laminar  to  turbu¬ 
lent  boundary-layer  transition  is  a  critical  parameter 
in  hypersonic  vehicle  design  [1],  The  transition  lo¬ 
cation  directly  affects  estimates  of  aeroheating  and 
skin  friction  drag,  which  in  turn  affect  heat  shield 
weight  and  materials,  range,  and  payload  capacity. 
Vehicle  designers  can  now  compute  the  laminar  and 
turbulent  aeroheating  with  good  accuracy;  in  many 
cases,  the  largest  uncertainty  in  calculating  the  total 
heat  flux  to  the  vehicle  results  from  the  estimates  of 
transition  location  [3]. 

Transition  Prediction  Methods 

Transition  prediction  for  vehicle  design  has  his¬ 
torically  relied  on  empirical  correlations  and  extrap¬ 
olations  from  wind  tunnel  experiments.  Many  corre¬ 
lations  exist,  and  most  fit  some  subset  of  the  exper¬ 
imental  data.  However,  no  empirical  model  can  ac¬ 
curately  predict  transition  for  a  general  dataset  [1]. 
In  fact,  Schneider  [4]  has  shown  that  for  a  general 
data  set  the  predictions  generated  by  different  cor¬ 
relations  scatter  by  a  factor  of  three  in  momentum 
thickness  Reynolds  number  and  an  order  of  magni¬ 
tude  in  arc  length  Reynolds  number.  Extrapolating 
from  wind  tunnel  experiments  is  also  problematic. 
Any  single  ground  facility  can  only  match  a  small 
subset  of  the  dozens  of  parameters  that  affect  tran¬ 
sition  [5].  Further  complicating  matters,  the  accous- 
tic  noise  produced  by  the  turbulent  boundary-layer 
on  the  walls  of  conventional  hypersonic  wind  tun¬ 
nels  produces  a  freestream  disturbance  environment 
markedly  different  from  that  of  flight.  These  distur¬ 
bances  have  been  shown  to  affect  not  only  the  tran¬ 
sition  location,  but  also  the  parametric  trends  [6]. 
This  has  led  to  a  push  for  the  development  of  ‘quiet’ 
tunnels,  but  there  are  currently  no  hypersonic  wind 
tunnels  that  are  quiet  at  significant  Reynolds  num¬ 
bers  anywhere  in  the  world  [7] .  Flight  tests  can  repli¬ 
cate  all  of  the  relevant  parameters,  but  only  limited 
information  is  available,  and  its  expense  is  normally 
prohibitive. 

Linear  Stability  Theory  (LST)  [8]  and  the 
Parabolized  Stability  Equations  (PSE)  [9]  have  been 
coupled  with  the  semi-empirical  method  to  pre¬ 
dict  hypersonic  transition  on  simple  geometries.  Re¬ 
cently,  Malik  [10]  showed  that  transition  locations 
could  be  correlated  on  two  hypersonic  flight  tests 
using  the  PSE  for  chemically  reacting  flows  and  the 

method.  Although  these  methods  have  shown 


success  in  research  settings  and  in  design  for  sub¬ 
sonic  and  transonic  applications,  hypersonic  vehicle 
designers  are  hesitant  to  stake  expensive  programs 
on  these  newer  methods  until  they  have  been  shown 
to  be  better  validated  and  simpler  to  use  [11]. 

Since  hypersonic  vehicles  remain  a  priority,  and 
current  design  methods  are  inadequate,  an  im¬ 
proved  method  for  hypersonic  boundary-layer  tran¬ 
sition  prediction  is  needed.  The  overall  project  that 
this  work  is  a  part  of  aims  to  provide  mechanism- 
based  methods  that  are  suitable  for  design  pur¬ 
poses,  including  predictions  using  the  PSE  and  the 

method.  This  will  not  be  the  final  solution  to 
the  problem  of  transition  prediction.  However,  it  is 
hoped  that  the  inclusion  of  more  of  the  relevant  flow 
physics  will  result  in  a  more  accurate  and  more  ro¬ 
bust  method  than  those  which  are  currently  in  use. 
The  present  work  is  concerned  primarily  with  the 
verification  and  validation  of  the  2D/axisymmetric 
version  of  PSE/e^  code  package. 

Verification  and  Validation 

Before  computational  results  can  be  used  to 
make  important  decisions,  measures  must  be  taken 
to  ensure  confidence  in  the  accuracy  of  the  predic¬ 
tions.  Although  this  issue  has  received  considerable 
attention  in  recent  years,  developments  in  this  area 
have  not  kept  pace  with  increases  in  code  size  and 
complexity  and  the  increasing  reliance  on  computer 
tions  to  reduce  the  number  and  scope  of  wind  tunnel 
and  flight  test  cases  [12]. 

Blottner  [13]  defines  code  validation  as  solving 
the  right  governing  equations  and  code  verification 
as  solving  the  governing  equations  right.  Examples 
of  validation  include  confirming  the  applicability  of 
the  assumptions  made  and  comparing  the  accuracy 
of  the  closure  models  throughout  the  relevant  par 
rameter  space.  Verification  means  evaluating  the 
accuracy  of  the  numerical  procedures  used  to  solve 
the  governing  equations.  Examples  of  verification  in¬ 
clude  checking  for  general  programming  errors  and 
ensuring  that  results  are  converged  and  grid  inde¬ 
pendent.  Code  validation  can  only  be  accomplished 
by  detailed  comparisons  to  experimental  data  [14], 
and  the  validation  of  applicable  equations  and  mod¬ 
els  is  a  responsibility  of  the  community  at  large.  In 
general,  a  code  must  be  verified  before  it  can  be  fully 
validated.  For  this  reason,  this  work  will  focus  pri¬ 
marily  on  verification,  although  the  validation  of  the 
transport  properties  will  be  discussed. 

Several  approaches  exist  for  code  verification,  in¬ 
cluding  grid  refinement  studies,  varying  the  bound¬ 
ary  conditions,  solving  problems  that  have  an  anar 
lytical  solution,  and  comparing  results  with  previ- 
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ously  verified  codes  [12,14].  Grid  refinement  studies 
and  comparisons  with  other  codes  will  be  used  in 
the  present  effort.  This  work  expands  on  previous 
verification  and  validation  by  Johnson  and  Johnson 
et  al.  [15-18].  Additional  details  about  the  work 
presented  here  are  in  Robarge  [19]. 

Code  Description 

The  computations  presented  in  this  paper  were  per¬ 
formed  using  the  STABL  code  package.  STABL  is 
designed  to  be  a  comprehensive  boundary-layer  sta¬ 
bility  analysis  and  transition  prediction  tool,  and  it 
includes  a  mean  flow  code,  a  stability  code,  and 
other  utilities  such  as  a  Graphical  User  Interface 
(GUI)  and  a  grid  generator.  STABL  is  modular  in 
nature,  allowing  individual  components  to  be  used 
instead  of  the  whole  package.  Further  details  about 
STABL  and  its  components  can  be  found  in  Johnson 
and  Candler  [18]  and  Robarge  [19]. 

The  laminar  mean  flows  were  obtained  using 
the  code  DPLR2D,  which  solves  the  axisymmet- 
ric  Navier-Stokes  equations  using  the  Data  Parallel 
Line  Relaxation  (DPLR)  method  [20]  with  modified 
Steger- Warming  flux  vector  splitting  [21].  The  solver 
is  second  order  accurate  in  both  the  streamwise  and 
wall-normal  directions,  although  it  drops  to  first  or¬ 
der  accuracy  near  the  shock.  DPLR2D  and  the  PSE- 
Chem  stability  code  can  both  handle  a  range  of  real 
gas  effects,  including  perfect  gas  or  reacting  flows 
with  chemical  non-equihbrium  using  a  general  set 
of  chemical  reactions.  Thermal  non-equilibrium  is 
treated  using  a  two-temperature  model  with  both 
translational  and  vibrational  modes.  Chemical  and 
thermal  equilibrium  can  both  be  simulated  through 
large  relaxation  rates. 

The  grids  used  for  the  mean  flow  solutions 
were  created  using  the  grid  generator  included  with 
STABL.  The  model  surface  and  outer  boundaries 
were  each  specified  using  a  single  arc  section  and  one 
or  more  line  segments.  In  the  case  of  a  sharp  cone, 
only  line  segments  were  used.  The  outer  boimdary 
was  specified  so  as  to  closely  match  the  local  shock 
angle.  The  lines  of  constant  i  were  constrained  to  be 
always  normal  to  tbe  wall.  The  surface  spacing  was 
uniform  on  the  nose,  and  exponential  stretching  us¬ 
ing  multiple  segments  distributed  points  smoothly 
on  the  cone  frusta.  The  lines  of  constant  j  were 
constructed  by  interpolating  between  spacings  spec¬ 
ified  at  the  upstream  and  downstream  boundaries. 
Exponential  stretching  was  used  to  cluster  points 
within  the  boundary-layer.  Experiments  with  clus¬ 
tering  points  at  the  wall  and  boundary-layer  edge 
showed  no  significant  effect  on  stability  results  [19]. 


The  code  PSE-Chem  is  used  to  analyze  the  sta¬ 
bility  of  the  laminar  mean-flow  profiles.  PSE-Chem 
uses  the  linear  Parabolized  Stability  Equations  cou¬ 
pled  with  the  method  to  predict  transition.  De¬ 
tails  of  the  PSE  can  be  found  in  Herbert  [9],  and 
the  full  deriviation  used  in  PSE-Chem  is  given  by 
Johnson  [17].  Linear  stability  theory  is  used  to  pro¬ 
vide  initial  wavenumber  guesses  for  the  PSE  march¬ 
ing  using  the  global  and  local  two-step  procedure  of 
Malik  [22]. 

The  disturbance  equations  axe  obtained  from 
the  Navier-Stokes  equations  by  first  decomposing 
the  instantcineous  flow  into  a  mean  and  a  fluctuat¬ 
ing  component,  substituting  into  the  Navier-Stokes 
equations,  and  subtracting  the  mean  flow  equations. 
The  disturbance  quantities  are  decomposed  into  a 
rapidly  varying  wave  function  and  a  slowly  varying 
shape  function,  as  shown  in  Equation  1. 

=  (1) 

The  shape  function  is  parabolized,  while  all  of  the 
ellipticity  is  retained  in  the  wave  function. 

Transition  prediction  is  accomplished  using  the 
method.  The  value  e^,  where  N  is  given  by 
Equation  2,  represents  the  total  growth  factor  of  a 
small  amplitude  initial  disturbance. 


In  linear  stability  theory,  only  the  contribution  of 
-Of  is  used  to  calculate  N  factors.  In  the  PSE,  the 
growth  rate  of  the  disturbances  includes  contribu¬ 
tions  from  both  the  imaginary  component  of  the 
streamwise  wavenumber  and  the  change  in  the  ki¬ 
netic  energy  of  the  shape  function,  as  shown  in  Equa¬ 
tions  3  and  4. 


E=  f  p{\uf  +  |t)p  +  |u)p)dC  (4) 
Jn 

Mack  [8]  showed  that  the  disturbances  can  have 
several  modes.  The  first  mode  is  analogous  to  the 
TS  waves  of  incompressible  flow.  The  first  mode  is 
damped  by  wall  cooling  and  is  most  amplified  when 
it  is  at  an  oblique  angle.  The  second  mode  can  be 
thought  of  as  a  trapped  acoustical  wave.  It  is  ampli¬ 
fied  by  wall  cooling,  and  it  is  most  amplified  when  it 
is  2D.  The  higher  modes  always  have  lower  amplifi¬ 
cation  than  the  first  or  second  mode,  and  they  were 
not  considered  here. 
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Table  1:  Test  conditions  for  Stetson  et  al.  [23]  blunt 
cone  experiment  for  all  computations  presented  here. 


Model  Specifications 

Nose  radius  (mm) 

3.81 

Half  angle 

7° 

Length  (m) 

1.016 

Measured  Flow  Conditions 

Moo 

7.99 

Po  (Pa) 

4.0x10® 

To  (K) 

750 

Pcx>  (Pa) 

410 

Inferred  Flow  Conditions 

Too  (K) 

54 

Poo  (kg/m^) 

0.0027 

Reoo/ft 

2.5x10® 

Second  Mode  Verification: 
Stetson’s  Blunt  Cone 

The  second  mode  calculations  in  STABL  were  veri¬ 
fied  by  benchmarking  agsdnst  other  codes  since  ac¬ 
curate  experimental  data  were  not  available  for  val¬ 
idation.  This  exercise  uncovered  a  number  of  issues, 
including  bugs,  alternative  algorithms,  and  different 
transport  property  models.  These  discoveries  led  to 
changes  that  improved  the  quality  of  STABL.  How¬ 
ever,  this  was  not  an  ideal  situation,  as  many  of  the 
details  of  the  other  calculations  are  unknown.  The 
methods  used  are  not  identical,  as  STABL  uses  the 
PSE  and  is  designed  for  chemically  reacting  fiows, 
whereas  all  of  the  benchmark  codes  use  LST  for  per¬ 
fect  gas  air  flows.  In  addition,  unknown  errors  may 
exist  in  the  other  codes  or  their  useage  that  could 
account  for  differences  in  the  results. 

Experimental  Conditions 

The  conditions  chosen  for  the  exercise  were 
those  of  the  blunt  cone  experiments  of  Stetson  et 
al.  [23]  with  the  3.81  mm  nose  radius.  Stetson  per¬ 
formed  hot  wire  measurements  in  the  boundary  layer 
of  a  sphere-cone  at  zero  degrees  angle  of  attack  and 
Mach  8  in  Tunnel  B  at  AEDC.  The  test  conditions 
for  all  Stetson  computations  in  this  paper  are  sum¬ 
marized  in  Table  1.  The  freestream  temperature 
and  pressure  are  required  inputs  for  STABL,  and 
they  were  calculated  from  the  specified  flow  condi¬ 
tions  using  the  isentropic  relations  and  the  perfect 
gas  law. 

The  Stetson  et  al.  experiment  has  been  analyzed 
by  numerous  other  researchers  using  many  different 


codes,  which  makes  it  an  ideal  benchmark  for  wind- 
tunnel  conditions.  Comparisons  were  made  with  re¬ 
sults  from  References  [24-32].  Most  of  these  cita¬ 
tions  are  summarized  by  Schneider  [3]. 

Schneider  [29]  has  shown  that  there  are  a  num¬ 
ber  of  open  issues  associated  with  this  experiment. 
The  axial  station  s/r„  =  175  is  used  for  the  vast 
majority  of  the  computations  presented  in  the  lit¬ 
erature.  However,  Stetson’s  data  show  considerable 
nonlinearity  at  this  station  which  will  not  be  cap¬ 
tured  by  a  LST  analysis.  Overshoots  are  present  in 
the  pitot  pressure  profiles,  which  are  thought  to  be 
due  to  interference  from  the  probe  and  have  not  ap¬ 
peared  in  any  computation.  Most  researchers  have 
assumed  the  wall  is  adiabatic,  but  the  measured  sur¬ 
face  temperature  is  about  20%  below  the  adiabatic 
temperature.  The  present  work  is  not  an  attempt  to 
resolve  the  discrepancies  between  computations  and 
experiments.  Rather,  the  intent  is  to  replicate  the 
computational  results  obtained  with  other  codes  to 
partially  verify  STABL.  The  large  body  of  indepen¬ 
dent  computational  data  makes  this  case  ideal  for 
verification,  but  these  unresolved  issues  make  com¬ 
plete  validation  from  the  experimental  data  impos¬ 
sible. 

A  number  of  issues  were  addressed  to  try  to  im¬ 
prove  the  agreement  between  the  stability  results 
of  STABL  and  those  of  other  researchers.  Discus¬ 
sion  of  the  laminar  transport  properties,  freestream 
conditions,  and  grid  generation  will  be  given  here. 
Robarge  [19]  describes  other  areas  that  were  consid¬ 
ered,  including  the  method  of  boundairy-layer  edge 
detection,  the  wall  temperature,  the  various  options 
available  in  STABL,  the  compilers  used,  variable 
normalization,  and  the  numerical  behavior  of  PSE- 
Chem.  These  issues  all  had  a  smaller  effect  on  the 
stability  results. 

Laminar  Transport  Properties 

The  lamineir  transport  properties  are  an  essen¬ 
tial  component  of  any  numerical  model.  Accepted 
models  exist  for  usage  within  normal  ranges  of  tem¬ 
perature  and  pressure.  However,  there  is  consider¬ 
able  subjectivity  in  determining  the  best  model  for 
very  high  or  low  temperatures,  and  the  best  choice 
for  one  extreme  is  generally  not  the  best  choice  for 
the  other.  Lyttle  and  Reed  [33]  showed  for  high  tem¬ 
peratures  that  the  transport  properties  employed 
can  have  a  large  impact  on  hypersonic  stability  re¬ 
sults. 

The  components  of  STABL  were  designed  pri¬ 
marily  for  high  enthalpy  air  flows,  such  as  those 
found  in  hypersonic  flight  or  in  shock  tunnels.  These 
flows  are  characterized  by  extremely  high  temper- 
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atures,  causing  dissociation,  ionization,  and  chemi¬ 
cal  and  thermal  non-equilibrium,  phenomena  known 
collectively  as  real  gas  effects.  Much  more  compli¬ 
cated  flow  models  are  required  to  accurately  model 
these  real  gas  effects  than  in  situations  where  perfect 
gas  behavior  can  be  assumed.  In  contrast,  the  wind 
tunnel  cases  discussed  here  all  have  edge  tempera¬ 
tures  below  100  K.  For  the  STABL  suite  to  be  valid 
in  the  low  enthalpy  regime,  transport  property  mod¬ 
els  appropriate  to  that  regime  must  be  employed. 

All  of  the  computations  cited  previously  were 
performed  with  a  Sutherland  viscosity  law.  Not  all 
of  the  references  provide  details  about  the  thermal 
conductivity  modeling,  but  those  that  did  specify 
the  method  assumed  a  perfect  gas  with  a  constant 
Prandtl  number.  It  is  probable  that  those  who  did 
not  specify  a  method  used  this  technique  as  well. 
Because  all  of  the  computations  in  the  literature 
treated  the  air  as  a  single  gas  rather  than  a  mix¬ 
ture,  none  included  the  effects  of  diffusion.  Due  to 
the  low  temperatures  involved,  the  production  of 
monatomic  species  is  negligible.  The  binary  diffu¬ 
sion  model  should  therefore  not  be  a  significant  fac¬ 
tor,  and  its  accuracy  was  not  investigated. 

Viscosity 

STABL  originally  used  Blottner’s  [34]  viscos¬ 
ity  model  for  reacting  flows  to  determine  individual 
species  viscosities  [17].  This  model  uses  empirical 
curve  flts  of  the  form  given  by  Equation  5  to  com¬ 
pute  individual  viscosities  in  kg/ra-s  for  each  species 
in  the  mixture. 

p.  =  (5) 

The  viscosity  of  the  mixture  is  found  by  Wilke’s 
semi-empirical  mixing  law. 

Two  sets  of  the  coefficients  A,,  B,,  and  C, 
are  given  in  Table  2.  The  high  temperature  set 
is  valid  over  the  temperature  range  1000  K-30000 
K  [34].  The  low  temperature  set  represents  unpub¬ 
lished  data  provided  by  Blottner  and  is  valid  over 
the  range  100  K-10000  K  [35].  Values  for  the  species 
NO,  N,  and  O  are  given  in  Reference  [34]  and  are  the 
same  for  both  temperature  ranges.  These  species 
have  mass  fractions  many  orders  of  magnitude  be¬ 
low  those  of  N2  and  O2  at  wind  tunnel  temperatures, 
and  their  viscosity  modeling  was  not  investigated. 

Figure  1  shows  a  comparison  of  the  two  meth¬ 
ods  given  in  Table  2  and  the  Sutherland  law  plotted 
with  the  experimental  data  of  Grieser  and  Goldth- 
waite  [36]  and  Matthews  et  al.  [37].  The  Blottner 
high  temperature  model  is  seen  to  differ  greatly  from 
the  experimental  data  at  low  temperatures.  At  80 


Table  2:  Blottner’s  Viscosity  Curve  Fit  Coefl&cients 


Coefficient 

High  Temp  [34] 

Low  Temp  [35] 

AjVa 

0.0268142 

0.0089993 

Bjvj 

0.3177838 

0.6039338 

Cjvj 

-11.3155513 

-12.4453814 

A02 

0.044929 

-0.0255541 

B02 

-0.0826158 

1.0503525 

C02 

-9.2019475 

-13.7080219 

K,  which  is  approximately  the  edge  temperature  for 
the  Stetson  case,  the  Blottner  model  overpredicts 
the  viscosity  by  70%.  The  Sutherland  law  diverges 
from  the  empirical  models  above  approximately  1900 
K,  as  is  expected  based  on  the  approximations  made 
in  its  derivation  [38]. 

Based  on  this  analysis,  STABL  was  modified 
to  incorporate  different  viscosity  models  for  each 
temperature  regime.  STABL  now  uses  the  Suther¬ 
land  law  for  low  temperatures,  the  high  tempera¬ 
ture  Blottner  law  for  high  temperatures,  and  the  low 
temperature  Blottner  law  for  the  intermediate  range. 
Blending  functions  are  used  to  ensure  a  smooth  tran¬ 
sition  between  the  laws.  The  resulting  method  is 
shown  as  the  STABL  Blended  curve  in  Figure  1. 
Agreement  with  the  experimental  data  is  consider¬ 
ably  better. 

It  is  important  to  note  that  none  of  the  three 
laws  implemented,  including  the  Sutherland  law,  are 
necessarily  accurate  below  100  K.  Several  models  ex¬ 
ist  for  this  low  temperature  regime.  The  Chapman- 
Enskog  model  [39]  is  based  on  the  kinetic  theory  of 
gases  and  employs  collision  integrals  to  model  the 
intermolecular  forces.  Keyes  [40]  and  Maitland  and 
Smith  [41]  both  proposed  empirical  curve  fits,  and 
Mack  [8]  proposed  a  model  hereafter  referred  to  as 
the  “Linear  Sutherland  Law”  that  is  a  linear  extrap¬ 
olation  of  the  Sutherland  Law  below  the  Sutherland 
constant.  Because  the  experimental  data  is  so  lim¬ 
ited  for  this  temperature  range,  there  is  no  clear-cut 
choice  for  the  best  model  to  use.  The  Sutherland  law 
was  used  for  all  computations  in  this  work  to  provide 
consistency  with  the  other  computations  referenced. 


Thermal  Conductivity 

STABL  calculates  thermal  conductivities  for 
both  translational  and  vibrational  modes  [17].  It 
uses  Eucken’s  relation,  given  in  Equations  6-8,  to 
calculate  the  translational  thermal  conductivity  of 
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Figure  1:  A  comparison  of  the  viscosity  computed 
with  both  sets  of  coefficients  of  the  Blottner  form, 
the  Sutherland  law,  and  experimental  data.  The 
line  labeled  STABL  Blended  represents  the  modi¬ 
fied  method  implemented  as  a  result  of  the  author’s 
analysis. 


Figure  2:  Thermal  conductivity  as  computed  by 
STABL  before  and  after  viscosity  law  modification 
and  with  a  constant  Prandtl  number  of  0.72  with 
constant  specific  heats  and  the  Sutherland  law.  Ex¬ 
perimental  data  are  shown  for  reference. 


each  species. 
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The  vibrational  thermal  conductivity  is  given  by 
Equation  9. 


kvib,  —  ftgCvfVib,  —  ftg 


de„. 

di; 


(9) 


The  mixture  thermal  conductivities  are  computed 
with  Wilke’s  semi-empirical  mixing  law. 

Figure  2  shows  a  comparison  of  the  thermal  con¬ 
ductivity  as  computed  by  STABL  and  as  computed 
using  a  constant  Prandtl  number  of  0.72.  The  curve 
marked  STABL  Original  uses  the  Blottner  high  tem¬ 
perature  viscosity  model,  the  curve  marked  STABL 
Blended  uses  the  blended  viscosity  model  described 
earlier,  and  the  constant  Prandtl  number  curve  uses 
a  Sutherland  Law.  Experimental  data  of  Taylor  and 
Johnston  [42],  Keyes  [43],  and  Vines  [44]  are  also 
presented.  The  STABL  blended  method  agrees  very 


well  at  low  temperatures  with  both  the  experimen¬ 
tal  data  and  the  constant  Prandtl  number  method. 
The  use  of  the  blended  viscosity  law  causes  the  error 
in  the  calculated  thermal  conductivity  with  respect 
to  the  experimental  data  to  drop  firom  70%  to  4% 
at  80  K.  However,  all  three  of  the  methods  differ 
significantly  from  the  experimental  data  at  higher 
temperatures.  This  may  be  due  to  the  importance  of 
the  vibrational  energy  mode  at  higher  temperatures, 
but  this  was  not  investigated  further.  Given  the 
close  agreement  between  the  constant  Prandtl  num¬ 
ber  method  and  that  used  by  STABL,  this  should 
not  be  a  source  of  difference  between  STABL  and  the 
computational  results  of  other  researchers.  However, 
it  could  have  an  impact  on  the  accuracy  of  actual 
stability  predictions. 

Figure  3  shows  the  effect  of  the  viscosity  model 
on  the  amplification  rate.  The  amplification  curve 
peak  is  shifted  from  (tu  kHz, -ai  l/m)=(136,13.4) 
to  (132,17.3),  which  is  an  increase  in  the  peak  am¬ 
plification  rate  of  29%  and  a  decrease  of  3%  in  the 
frequency  at  which  that  occurs.  The  decrease  in  the 
frequency  causes  better  agreement  with  the  other  re¬ 
searchers,  but  the  increase  in  the  amplification  rate 
causes  worse  agreement.  Although  this  improves 
agreement  in  the  mean  flow  profiles,  the  agreement 
in  the  ampliflcation  curve  is  not  improved  by  the 
change  in  the  viscosity  model. 
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Figure  3;  Amplification  rate  with  the  original  and 
blended  viscosity  models  for  Stetson’s  blunt  cone  at 
s/r„  =  175  with  r„  =  3.81  mm  using  the  conditions 
given  in  Table  1. 


Table  3:  Unit  Reynolds  number  for  the  Stetson  et 
al.  blunt  cone  experiment  calculated  using  various 
viscosity  models. 


Viscosity  Model 

RCoo/ft 

Sutherland 

Linear  Sutherland 
Keyes 

Chapman-Enskog 

3.546  X  10“® 
3.771  X  10"® 
3.893  X  10-® 
3.788  X  10"® 

2.69  X  10® 
2.53  X  10® 
2.45  X  10® 
2.53  X  10® 

Input  Conditions 

In  Reference  [23],  Stetson  et  al.  specify  the  quan¬ 
tities  Moo,  Po,  Toi  Poo,  and  unit  Reynolds  number. 
The  Mach  number,  pressures,  and  stagnation  tem¬ 
perature  are  measured  quantities.  The  freestream 
temperature,  which  is  a  required  input  for  many 
codes,  including  STABL,  can  be  calculated  using 
the  isentropic  relations.  For  this  case,  it  is  ap¬ 
proximately  54.4  K.  The  unit  Reynolds  number  is 
a  derived  quantity  that  depends  on  the  viscosity 
model,  which  can  be  a  source  of  uncertainty  at  that 
freestream  temperature.  Table  3  shows  a  compar¬ 
ison  of  the  unit  Reynolds  number  calculated  using 
different  viscosity  laws  in  common  use  for  the  low 
temperature  regime. 

It  seems  likely  that  Stetson  et  al.  used  either 
the  linear  Sutherland  law  or  the  Chapman-Enskog 


Figure  4;  Amplification  rates  computed  using  two 
different  values  of  freestream  density  for  Stetson’s 
blunt  cone  at  s/r„  =  175  with  r„  =  3.81  mm  using 
the  conditions  given  in  Table  1. 


model  to  calculate  the  unit  Reynolds  number  they 
reported.  However,  the  other  researchers  employed 
a  standard  Sutherland  law  for  their  computations, 
and  the  unit  Reynolds  number  is  an  input  condi¬ 
tion  for  many  CFD  codes.  To  truly  match  the 
freestream  conditions  of  Stetson’s  experiment  using 
a  Sutherland  law  code,  a  unit  Reynolds  number  of 
2.69  X  10®/ft  would  need  to  be  used,  rather  than  the 
2.50  X  10®/ft  specified  in  Reference  [23]. 

The  STABL  mean  flow  solver  uses  the 
freestream  temperature,  freestream  Mach  number, 
and  freestream  density  ^  input  conditions.  The 
freestream  density  is  0.02652^  when  computed 
from  the  stagnation  temperature  and  pressure  us¬ 
ing  the  perfect  gas  law  and  the  isentropic  relations. 
If  a  Sutherland  viscosity  law  is  used  with  a  unit 
Reynolds  number  of  2.50  x  10®/ft,  the  freestream 
density  becomes  0.02464^.  The  first  value  appears 
to  be  more  physically  accurate,  and  should  be  used 
for  comparison  with  the  experimental  data.  How¬ 
ever,  the  second  number  is  most  likely  what  was  ac¬ 
tually  used  in  the  other  researchers’  computations, 
and  thus  it  should  be  used  for  code  verification  with 
the  current  benchmarks.  Figure  4  shows  the  ef¬ 
fect  of  changing  the  density  input.  Use  of  the  unit 
Reynolds  number  based  value  shifts  the  amplifica¬ 
tion  rate  curve  closer  to  the  other  researchers’  re¬ 
sults,  moving  the  peak  from  (132,17.9)  to  (127,16.6). 
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Figure  5:  N  factors  computed  using  mean  flows  with 
varying  grid  resolutions  for  Stetson’s  blunt  cone  with 
r„  =  3.81  mm  using  the  conditions  given  in  Table  1. 


Grid  Generation 

The  requirement  for  accurate  second  derivatives 
of  the  velocity  and  temperature  profiles  through¬ 
out  the  boundary-layer  means  that  much  finer  grids 
must  be  used  for  stability  analysis  than  for  most 
other  CFD  applications.  Figure  5  shows  the  N  fac¬ 
tors  computed  using  several  grids  with  varying  res¬ 
olution.  The  five  curves  shown  used  grids  with  106, 
121,  140,  239,  and  276  points  in  the  boimdary-layer. 
The  grids  marked  (b)  in  Figure  5  had  an  outer 
boundary  that  better  matched  the  shock  angle  than 
those  marked  (a).  The  grids  marked  (b)  produce 
slightly  smaller  N  factors  than  the  grids  marked  (a). 
There  is  no  discernable  trend  with  the  total  grid  size. 
This  suggests  that  the  differences  are  due  to  the  sen¬ 
sitivity  of  the  stability  calculations  to  small  grid  vari¬ 
ations,  and  it  cannot  be  conclusively  said  that  one 
grid  shape  is  better  than  the  other.  The  scatter  of 
approximately  5%  at  s  =  1.0  m  may  represent  an 
uncertainty  in  the  residts  that  cannot  be  eliminated 
through  increasing  mesh  density.  In  any  case,  the 
differences  are  slight  compared  with  the  overall  im- 
certainty  associated  with  the  use  of  the  method, 
and  sufficient  mesh  spacing  is  used  for  the  present 
results. 

Results 

The  mesh  used  for  the  final  computation  con¬ 
tained  440  axial  points  and  350  body-normal  points. 


Figure  6:  Mean  velocity  profile  for  conditions  of 
Stetson  et  al.  [23]  at  s/r„  =  175  with  r„  =  3.81 
mm  compared  with  the  results  of  Kufner  et  al.  [26]. 


Figure  7:  Second  normal  derivative  of  velocity  profile 
for  conditions  of  Stetson  et  al.  [23]  at  s/r„  =  175 
with  r„  —  3.81  mm  compared  with  the  results  of 
Ku&ier  et  al.  [26]  and  Esfahanian  [25]. 
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Figure  8:  Amplification  rate  for  conditions  of  Stet¬ 
son  et  al.  [23]  at  s/r„  =  175  with  Vn  =  3.81  mm. 


and  it  corresponds  to  the  440x350  (b)  line  in  Fig¬ 
ure  5.  Chemical  and  thermal  nonequilibrium  effects 
were  modeled,  and  a  freestream  air  mixture  with 
standard  mass  fractions  of  76.7%  N2  and  23.3%  O2 
was  assumed.  The  freestream  density  was  set  to 
0.02464^  to  match  the  conditions  of  the  other  re¬ 
searchers.  Thermal  equilibrium  was  assumed  in  the 
freestream.  STABL  version  1.26  was  used  for  both 
the  mean  flow  and  the  stability  calculations. 

Figure  6  shows  the  mean  tangential  flow  veloc¬ 
ity  at  s/r„  =  175.  Data  in  Figure  la  of  Kufher 
et  al.  [26]  was  scanned  and  digitized  and  is  plotted 
on  Figure  6  for  comparison.  Agreement  is  gener¬ 
ally  good,  although  there  is  a  small  difference  near 
the  edge  of  the  boundary-layer.  Figure  7  shows  the 
second  derivative  of  the  velocity  profile  with  respect 
to  the  body-normal  distance  at  the  same  station. 
Data  from  Figure  lb  in  Reference  [26]  is  also  plotted. 
Again,  agreement  is  good  except  near  the  boundary- 
layer  edge.  Because  the  instabilities  arise  near  the 
boundary-layer  edge,  these  differences  could  be  sig¬ 
nificant. 

Figure  8  shows  the  amplification  curve  at  s/r„  = 
175  compared  with  that  computed  by  other  re¬ 
searchers.  The  peak  frequency  is  slightly  lower  than 
the  mean  peak  frequency,  but  it  is  within  the  scat¬ 
ter  of  the  other  results.  The  peak  amplification  rate 
is  higher  than  any  of  the  other  researchers.  The 
magnitude  of  the  difference  in  the  peak  amplifica¬ 
tion  rate  is  comparable  to  the  scatter  in  the  other 


Table  4:  Location  of  the  peak  amplification  as  com¬ 
puted  by  STABL  and  several  researchers. 


Computation 

w  (kHz) 

-Oi  (1/m) 

STABL 

127 

16.5 

Malik  [24] 

132 

13.8 

Esfahanian  [25] 

129 

15.1 

Kufner  [26] 

126 

13.5 

Stilla  [26] 

122 

14.3 

Lyttle  [31] 

133 

14.7 

Zhong  [32] 

134 

15.2 

Figure  9:  Maximum  N  factors  obtained  for  various 
frequencies  for  the  Stetson  et  al.  blunt  cone  case  with 
r„  =  3.81  mm. 


results  shown.  Table  4  summarizes  the  location  of 
the  peaks  in  Figure  8.  A  comparison  with  the  results 
of  Malik  et  al.  [24]  at  s/r„  =  215  showed  the  same 
trends. 

Figure  9  displays  the  maximum  N  factors  ob¬ 
tained  for  various  frequencies  for  the  Stetson  et 
al.  blunt  cone  case.  To  obtain  these  N  factors,  a 
test  matrix  was  constructed  consisting  of  all  combi¬ 
nations  of  the  starting  locations  s  =  0.1, 0.15, . . . ,  1.0 
and  frequencies  ui  =  50, 60, ... ,  200  kHz.  PSE 
marching  was  conducted  at  each  point  in  the  test 
matrix,  and  the  maximum  N  factors  for  each  fre¬ 
quency  are  plotted  in  Figure  9.  K,  following  the 
work  of  Malik  et  al.  [24],  Stilla  [27],  and  Esfaha- 
nian  and  Hejranfar  [30],  transition  is  assumed  to  oc¬ 
cur  at  N  =5.5  in  this  noisy  environment,  STABL 
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Table  5:  Transition  location  predicted  by  STABL  for 
Stetson  et  al.  blunt  cone  case  compared  with  that 
predicted  by  other  researchers.  Transition  was  as¬ 
sumed  at  N=5.5  for  all  computations. 


Computation 

Str  (m) 

STABL 

1.00 

8.192  X  10® 

Esfahanian  PNS  [30] 

1.00 

8.21  X  10« 

Esfahanian  IPNS  [30] 

0.973 

7.98  X  10« 

Stilla  [27] 

0.956 

7.84  X  10« 

Malik  et  al.  [24] 

0.956 

7.84  X  10® 

5:6 — 
s(m) 


Figure  10:  Comparison  of  N  factors  calculated  by 
STABL  with  local  N  factors  of  Rosenboom  et  aJ.  [28] 
for  the  Stetson  et  al.  [23]  case  using  the  conditions 
given  in  Table  1. 


would  predict  transition  at  s  =  1.00  m  with  130 
kHz  the  most  amplified  frequency.  Table  5  shows  a 
comparison  of  the  transition  location  predicted  by 
STABL  with  that  of  the  other  researchers.  The  unit 
Reynolds  number  specified  by  Stetson  et  al.  of  2.50 x 
10®/ft  was  used  to  convert  between  the  and 

Str-  Agreement  with  the  other  researchers  is  better 
than  5%. 

Figure  10  shows  a  comparison  of  the  N  factors 
computed  by  STABL  with  the  local  N  factors  com¬ 
puted  by  Rosenboom  et  al.  [28].  Agreement  in  the 
location  of  amplification  of  the  various  frequencies  is 
fair,  with  no  clear  trends  evident.  The  STABL  140 
kHz  and  160  kHz  line  show  amplification  through¬ 
out  the  full  surface  length,  in  contrast  to  the  Rosen¬ 
boom  results,  which  begin  to  decay  at  s  =  0.85  and 


8  =  0.55  m,  respectively.  The  reason  for  the  dif¬ 
ferences  is  not  clear.  The  Rosenboom  et  al.  results 
are  not  given  in  Table  5  because  the  data  do  not  ex¬ 
tend  beyond  s  =  1.00  m  or  =  5.0,  but  it  is  clear 
that  the  predicted  transition  location  would  be  aft 
of  s  =  1.00  m. 

First  Mode  Verification 
Sharp  Cone  at  Mach  3.5 

A  sharp  cone  was  analyzed  using  the  conditions 
shown  in  Table  6  to  verify  STABL’s  ability  to  cal¬ 
culate  first  mode  disturbances.  The  conditions  cor¬ 
respond  to  experimental  data  obtained  by  Beckwith 
et  al.  [45]  in  the  NASA  Langley  Mach  3.5  pilot  quiet 
tunnel.  This  tunnel  is  designed  to  give  freestream 
disturbance  levels  comparable  to  flight.  The  condi¬ 
tions  chosen  correspond  to  the  fourth  unit  Reynolds 
number  data  set  in  Run  5  of  Reference  [45].  The 
first  mode  is  expected  to  be  dominant  at  this  low 
freestream  Mach  number  with  an  adiabatic  wall. 

This  case  was  previously  analyzed  by  Malik  as 
case  QTl  in  Reference  [46]  using  the  COSAL  sta¬ 
bility  code.  This  older  LST-based  code  calculates 
temporal  stability  and  uses  Caster’s  group  veloc¬ 
ity  transformation  to  obtain  spatial  stability  results. 
The  me£in  flow  was  obtained  from  a  boundary-layer 
code  and  provided  as  an  input  to  COSAL. 

Roger  Kinunel  [47]  performed  additional  compu¬ 
tations  for  this  case  using  the  stability  solver. 

This  more  recent  code  is  also  widely  used  for  hy¬ 
personic  stability  euaalysis.  Spatial  stabifity  results 
are  calculated  directly,  and  a  perfect  gas  model  is 
assumed.  For  the  cases  shown,  the  mean  flow  was 
provided  by  a  similarity  solver  built  into  the  e*^®***’ 
code.  These  computations  were  performed  at  the 
request  of  the  author,  and  a  good  deal  of  collabo¬ 
ration  and  discussion  with  Dr.  Kimmel  occurred  to 
troubleshoot  differences  between  the  results. 

A  grid  with  450  axial  points  and  350  normal 
points  was  used  for  the  DPLR2D  mean  flow  solu¬ 
tion.  Exponentied  stretching  was  used  in  both  direc¬ 
tions  to  cluster  grid  points  near  the  cone  tip.  This 
clustering  increased  the  number  of  points  within  the 
boundary-layer  and  reduced  the  degree  of  numeri¬ 
cal  error  caused  by  the  discontinuity  at  the  nose. 
The  number  of  grid  points  within  the  boundary- 
layer  varied  from  148  near  the  nose  tip  to  109  at 
the  base.  Chemical  and  vibrational  nonequilibrium 
effects  were  modeled  with  thermal  equilibrium  as¬ 
sumed  in  the  freestream,  and  the  blended  viscosity 
model  was  used  for  all  calculations. 
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Table  6:  Test  conditions  for  sharp  cone  at  Mach  3.5 


Condition 

Value 

Cone  half  angle 

5° 

Cone  length  (m) 

0.381 

Wall  temperature 

Adiabatic 

Moo 

3.5 

Po  (kPa) 

525 

To  (K) 

319.0 

Too  (K) 

92.53 

Poo  (kg/m^) 

0.02592 

Reoo/m 

2.74x10^ 

Sir  (m) 

0.278 

Figure  11  shows  wall-tangential  mean  flow  veloc¬ 
ity  profiles  computed  by  STABL  and  at  three 
axial  locations.  The  profiles  computed  by  STABL 
are  slightly  fuller  than  those  computed  by  for 

all  distances.  Figure  12  shows  the  mean  temperature 
profiles  for  the  same  locations.  Small  but  significant 
differences  are  seen  at  all  three  locations.  These  vis¬ 
ible  differences  in  the  temperature  and  velocity  can 
be  expected  to  have  a  significant  effect  on  the  sta¬ 
bility  results.  Meaningful  second  derivative  compar- 
isions  were  not  available  due  to  the  lack  of  precision 
in  the  output  of  the  e^®***'  similarity  solution. 

Several  mean  flow  solutions  using  different  num¬ 
bers  of  grid  points  within  the  boundary-layer  were 
analyzed  to  ensure  grid-independence  of  the  solu¬ 
tion.  Figure  13  shows  the  effect  of  the  different  grids 
on  the  N  factor  calculations.  These  calculations  are 
for  a  mode  at  w  =  104  kHz  and  P  =  2400/m  with 
marching  beginning  at  s  =  0.0339  m.  Four  of  the  so¬ 
lutions  are  very  similiar;  the  only  difference  is  seen 
for  the  grid  with  72  points  in  the  boundary-layer. 
The  reason  for  this  repeatable  difference  is  not  clear. 
The  grid  with  109  points  in  the  boundary-layer  was 
used  for  all  other  calculations  in  this  section. 

To  obtain  the  maximum  N  factors  for  this  case, 
a  test  matrix  was  constructed  consisting  of  w  = 
25, 50,..., 250  kHz,  P  =  1000, 2000, ..., 8000  1/m, 
and  starting  location  s  =  0.02, 0.04,  ...,0.12  m.  All 
of  the  480  combinations  of  these  three  parameters 
were  analyzed  in  approximately  16  hours  on  two  2.4 
GHz  processors.  The  modes  producing  the  maxi¬ 
mum  N  factors  at  each  location  are  shown  in  Fig¬ 
ure  14.  Amplification  began  at  s  =  0.02  m  for  the 
higher  frequencies  and  shortly  thereafter  for  the  75 
kHz  mode.  At  the  experimental  transition  location 
of  s  =  0.278  m,  the  most  amplified  mode  has  a  fre¬ 


Figure  11:  Mean  flow  velocity  profiles  for  the  sharp 
cone  at  Mach  3.5  at  several  axial  locations. 


quency  of  75  kHz,  a  ^  of  2000  1/m,  and  a  N  factor 
of  12.2.  This  N  factor  is  larger  than  the  r^lnge  of 
9-11  commonly  calculated  at  experimentally  mea¬ 
sured  transition  locations  in  low  disturbance  envi¬ 
ronments. 

Comparisons  were  made  between  the  present  re¬ 
sults  and  those  obtained  by  Malik  using  COSAL  [46] 
and  Kimmel  using  [47].  Malik  calculated  N 

factor  growth  rates  for  six  frequencies  ranging  from 
22-188  kHz.  His  most  amplified  mode  at  the  ex¬ 
perimental  transition  location  has  a  frequency  of  78 
kHz,  a  wave  angle  of  approximately  60°,  and  a  N 
factor  of  10.1.  Compared  to  the  transition  N  factor 
calculated  by  STABL  of  12.2,  this  is  a  significant 
difference,  although  the  frequency  and  wave  angle 
of  the  most  amplified  disturbances  agree  well. 

Figure  15  shows  a  comparison  of  the  N  factors 
calculated  by  STABL  and  those  calculated  by  Kim¬ 
mel  [47]  using  the  e^®'**  code  for  a  mode  at  104 
kHz  with  P  =  2400  1/m.  The  N  factors  calculated 
by  STABL  are  considerably  higher  than  those  calcu¬ 
lated  by  e^®***.  The  agreement  is  much  worse  than 
that  seen  in  Figures  11  and  12  for  the  mean  flow 
profiles.  The  difference  in  N  factor  may  be  due  to 
the  subtle  differences  in  the  mean  flow,  or  it  may 
be  due  to  differences  in  the  stability  solvers.  Plans 
for  future  work  include  using  the  mean  flow  gener¬ 
ated  by  the  similarity  solver  in  e^®***  for  stability 
calculations  by  PSE-Chem. 
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Figure  12:  Mean  flow  temperature  profiles  for  the 
sharp  cone  at  Mach  3.5  at  several  axial  locations. 


Figure  13:  N  factor  calculations  for  the  sharp  cone  at 
Mach  3.5  using  five  separate  mean  flow  grids.  Cal¬ 
culations  are  for  a  mode  with  w  =  104  kHz  and 
P  =  2400/m,  beginning  at  s  =  0.0339  m. 


Figure  14:  Maximum  N  factors  obtained  from  many 
combinations  of  frequency,  spanwise  wavenumber, 
and  starting  location  for  the  sharp  cone  at  Mach 
3.5. 


Figure  15:  N  factors  calculated  using  STABL  com¬ 
pared  with  results  of  Kimmel  [47]  using  for  a 

mode  at  104  kHz  with  P  =  2400  1/m  for  the  sharp 
cone  at  Mach  3.5. 
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Table  7:  Test  conditions  for  cone  with  r„  =  42.67 
mm  at  Mach  8 


Condition 

Value 

Nose  radius  (mm) 

42.67  mm 

Cone  half  angle 

7° 

Maximum  arc  length  (m) 

11.35 

Wall  temperature 

Adiabatic 

Moo 

7.99 

Po  (kPa) 

4000 

Tc,(K) 

750 

Too  (K) 

54 

Poo  (kg/m®) 

2.7x10-2 

Rcoo/ni 

2.5  xlO® 

High  Bluntness  Cone  at  Mach  8 

For  a  blunt  cone,  when  the  nose  radius  increases 
for  a  fixed,  high  freestream  Mach  number,  the  first 
mode  may  become  dominant.  For  this  reason,  a  very 
blunt  cone  was  analyzed  using  the  conditions  shown 
in  Table  7  to  further  verify  first  mode  computations. 
The  flow  conditions  match  those  of  the  Stetson  et 
al.  [23]  blunt  cone  experiment,  but  the  bluntness  was 
increased  while  the  ratio  of  the  nose-tip  radius  to  the 
body  length  was  kept  constant,  making  this  case  im¬ 
practical  for  wind  tunnel  experiments.  Rosenboom 
et  al.  [28,48]  analyzed  this  case  and  found  first  mode 
N  factor  growth  at  a  level  comparable  to  second 
mode  growth. 

An  accurate  mean  flow  was  computed  using  a 
mesh  with  450  axial  points  and  350  wall-normal 
points.  Approximately  100  grid  points  were  con¬ 
tained  within  the  boundary-layer  over  the  cone  frus¬ 
tum.  The  edge  Mach  number  was  approximately 
4  near  the  cone  shoulder  and  increased  to  approxi¬ 
mately  6.3  far  downstream  of  the  nose. 

Computations  were  performed  with  STABL  to 
compare  with  the  results  of  Rosenboom  et  al.  Fig¬ 
ure  16  shows  samples  of  the  first  mode  N  fac¬ 
tors  calculated  by  STABL.  A  test  matrix  consist¬ 
ing  of  combinations  of  s  =  0.5, 1.0, . . . ,  10.5  and 
^  =  20, 40, . . . ,  100  1/m  was  analyzed.  Combina¬ 
tions  involving  additional  values  of  ^  <  20  1/m  were 
also  analyzed  separately,  but  no  significant  amplifi¬ 
cation  was  found.  For  all  cases  shown,  a  frequency 
range  of  1-10  kHz  was  analyzed  at  the  starting  loca¬ 
tion  using  LST,  and  the  lowest  unstable  frequency 
was  used  for  PSE  marching.  For  these  conditions, 
the  marching  frequency  was  always  the  lowest  for 
which  converged  solutions  could  be  obtained.  P 
and  w  were  fixed  as  the  marching  progressed  down¬ 


Figure  16:  First  mode  N  factors  for  the  large  blunt¬ 
ness  cone  at  Mach  8. 


stream.  N=  5.5  is  first  reached  at  s  =  9.7  m  by  a 
3.33  kHz  instability  at  70°.  Despite  differences  in  the 
frequency  and  wave  angle,  the  transition  locations 
predicted  by  STABL  and  Rosenboom  et  al.  agree 
extremely  well.  In  addition,  the  general  characteris¬ 
tics  of  the  N  factor  curves  calculated  by  STABL  and 
Rosenboom  et  al.  agree  well.  Both  show  an  initial 
hump,  followed  by  brief  damping  and  a  long  rise. 
Table  8  summarizes  the  comparison. 

Figure  17  shows  the  second  mode  N  factor 
curves  generated  by  STABL.  A  frequency  range  of 
20-60  kHz  was  analyzed,  but  no  eigenvalues  were 
found  by  the  global  solver  at  frequencies  greater 
than  35  kHz  for  any  location.  STABL  selected  the 
critical  frequency  at  each  starting  location  to  begin 
PSE  marching.  Amplification  begins  at  s  =  6.5  m, 
slightly  behind  the  critical  location  of  s  =  5.5  m  cal¬ 
culated  by  Rosenboom  et  al.  An  N  factor  of  5.5  is 
first  reached  at  s  =  9.0  m  by  a  34  kHz  wave.  This  is 
slightly  ahead  of  the  location  calculated  by  Rosen¬ 
boom  of  s  =  9.6  m  for  the  local  calculation  or  s  =  9.8 
for  the  nonlocal  calculation,  but  the  difference  is  on 
the  order  of  the  starting  distance  resolution.  Overall 
second  mode  agreement  between  STABL  and  Rosen¬ 
boom  is  good.  Table  8  summarizes  the  comparison. 

An  outstanding  question  in  stability  theory  is 
the  nature  of  the  interaction  between  the  first  and 
second  modes.  Figure  18  shows  the  N  factors  for  the 
most  amplified  first  and  second  mode  waves  together 
with  the  edge  Mach  number  distribution.  The  dis¬ 
continuity  in  the  edge  Mach  number  at  s  =  3  m  is 
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Figure  17:  Second  mode  N  factors  for  the  large 
bluntness  cone  at  Mach  8. 


Table  8:  Transition  location  predictions  for  the  large 
bluntness  cone.  Transition  is  assumed  at  IV  =  5.5. 


Computation 

Sir  (m) 

u  (kHz) 

V' 

First  Mode 

STABL 

Km 

3.33 

Rosenboom  [28,48] 

2 

Second  Mode 

STABL 

9.0 

MM 

Rosenboom  (local) 

9.6 

msM 

IS 

Rosenboom  (nonlocal) 

9.8 

34 

IS 

caused  by  numerical  difficulties  in  defining  the  edge, 
as  discussed  further  in  Robarge  [19].  The  first  mode 
amplifies  over  the  full  length  of  the  body.  The  sec¬ 
ond  mode  begins  to  amplify  when  the  edge  Mach 
number  is  5.7,  and  it  becomes  dominant  when  the 
edge  Mach  number  is  6.0. 

The  calculations  also  showed  another  unstable 
mode  that  was  not  expected.  Figure  19  shows  the 
N  factors  calculated  by  two-dimensional  waves  at  1 
kHz.  The  different  curves  represent  different  starting 
locations.  The  waves  amplify  very  slightly  until  s  = 
3  m,  at  which  point  they  amplify  very  rapidly  until 
s  =  4  m.  They  then  decay  slightly  until  s  =  10 
m,  at  which  point  they  again  amplify  rapidly  until 
the  end  of  the  body.  N  =  5.5  is  first  reached  at 
8  =  10.6  m.  These  waves  were  found  serendipitously 


Figure  18:  A  comparison  of  the  largest  first  and  sec¬ 
ond  mode  N  factors  calculated  by  STABL  for  the 
large  bluntness  cone  at  M£ich  8. 


Figure  19:  N  factors  of  a  two-dimensional  instability 
at  1  kHz.  The  four  lines  represent  different  begin¬ 
ning  marching  locations. 
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during  a  search  over  the  1-10  kHz  frequency  range 
when  the  lowest  amplified  frequency  was  chosen  for 
marching.  Although  this  search  was  run  every  0.5 
m,  only  the  starting  locations  between  1.5  and  3.0  m 
exhibited  this  behavior.  The  reason  for  the  different 
amplification  between  the  different  starting  locations 
is  unknown.  This  might  be  a  new  phenomenon,  but 
further  investigation  is  necessary  before  it  can  be 
characterized. 

Conclusions 

Improved  methods  of  hypersonic  boundary-layer 
transition  prediction  require  extensive  verification 
and  validation.  Computations  have  been  conducted 
using  STABL  for  several  benchmark  cases. 

The  Stetson  et  al.  [23]  experiment  is  the  best- 
known  example  of  measured  second  mode  wave 
growth  on  a  blunt  cone.  This  effort  exposed  the 
inadequacy  of  the  method  originally  implemented  in 
STABL  for  computing  the  viscosity  of  low  tempera¬ 
ture  flows,  such  as  those  commonly  found  in  hyper¬ 
sonic  wind  tunnels.  Modifications  that  were  made  to 
use  a  blended  viscosity  law  helped  to  extend  the  ap¬ 
plicability  of  STABL  to  the  low  temperature  regime. 
A  related  issue  was  the  use  of  the  proper  freestream 
density.  Many  CFD  codes  require  a  xmit  Reynolds 
number  input,  but  the  unit  Reynolds  number  re¬ 
ported  in  Stetson  et  al.  is  not  consistent  with  the 
other  reported  quantities  if  a  Sutherland  viscosity 
law  is  used.  Thus,  any  researcher  that  used  the  unit 
Reynolds  number  specified  by  Stetson  et  al.  as  an  in¬ 
put  to  a  code  that  uses  the  Sutherland  law  modeled 
the  flow  conditions  of  the  experiment  incorrectly. 

Although  the  transition  location  predicted  by 
STABL  agreed  well  with  that  predicted  by  other  re¬ 
searchers,  there  is  considerable  scatter  in  the  ampli¬ 
fication  curves  Ccdculated  for  the  3.81  mm  Stetson 
cone  at  s/r„  =  175.  STABL  calculates  a  peak  am¬ 
plification  rate  that  is  more  unstable  them  any  of 
the  data  sets  used  for  comparison.  However,  the 
difference  is  approximately  equal  to  the  scatter  in 
the  other  results.  Furthermore,  when  the  viscosity 
model  and  method  of  calculating  the  freestream  den¬ 
sity  were  changed,  the  results  calculated  by  STABL 
shifted  by  amounts  similar  to  the  overall  scatter  in 
the  results.  This  shows  a  high  sensitivity  to  the 
physical  models  employed,  and  work  should  continue 
in  validating  these  models  for  various  flow  regimes. 

Computations  were  also  conducted  to  verify  the 
ability  of  STABL  to  predict  first  mode  disturbances. 
A  sharp  cone  at  Mach  3.5  and  a  very  blunt  cone 
at  Mach  8  were  analyzed.  The  sharp  cone  compar¬ 


ison  showed  good  agreement  in  the  transition  loca¬ 
tion  with  the  results  of  Malik  [46].  However,  the  N 
factor  computed  for  a  single  mode  differed  by  30% 
compared  to  a  new  corresponding  calculation  by  the 
^Mahk  code.  Future  plans  call  for  the  use  of  the 
mean  flow  solution  produced  by  in  PSE-Chem 
stability  calculations.  This  will  help  to  localize  the 
cause  of  the  difference.  Analysis  of  a  large-bluntness 
cone  showed  first  mode  amplification  over  most  of 
the  body.  Second  mode  amplification  began  farther 
aft,  but  it  soon  had  greater  total  amplitude  than 
the  first  mode.  The  predicted  transition  locations 
agreed  well  with  the  results  of  Rosenboom  et  al.  [28] 
for  both  the  first  and  second  mode. 

A  great  deal  of  progress  has  been  made  in  trans¬ 
forming  STABL  from  a  research  code  to  a  design 
tool.  The  computations  presented  here  help  to  pro¬ 
vide  confidence  in  the  ability  of  STABL  to  ana¬ 
lyze  boundary-layer  stability  and  predict  transition. 
These  results  can  also  serve  as  a  benchmark  to  al¬ 
low  other  researchers  using  STABL  to  verify  their 
method  of  operation.  Additional  computations  on 
different  cases  in  a  variety  of  flow  regimes  should  he 
performed  to  further  this  effort. 
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